Study of Intrinsic Spin Hall Effect and Orbital Hall Effect 
in Ad- and 5d- Transition Metals 



00 
O 

o 

(N 

X> 
<D 
PL, 

\o 

(N 



•4— > 

a 

c 

o 

> 
m 



O 



T. Tanaka 1 , H. Kontani 1 , M. Naito 1 , T. Naito 1 , D.S. Hirashima 1 , K. Yamada 2 and J. Inoue 3 
1 Department of Physics, Nagoya University, Furo-cho, Nagoya 464-8602, Japan. 
2 Engineering, Ritsumeikan University, 1-1-1 Noji Higashi, Kusastu, Shiga 525-8577, Japan. 
s Department of Applied Physics, Nagoya University, Furo-cho, Nagoya 464-8602, Japan. 

(Dated: February 26, 2008) 

We study the intrinsic spin Hall conductivity (SHC) in various 5<i-transition metals (Ta, W, Re, 
Os, Ir, Pt, and Au) and 4d-transition metals (Nb, Mo, Tc, Ru, Rh, Pd, and Ag) based on the 
Naval Research Laboratory tight-binding model, which enables us to perform quantitatively reliable 
analysis. In each metal, the obtained intrinsic SHC is independent of resistivity in the low resistive 
regime (p < 50^f2cm) whereas it decreases in proportion to p~ 2 in the high resistive regime. In 
the low resistive regime, the SHC takes a large positive value in Pt and Pd, both of which have 
approximately nine d-electrons per ion {rid = 9). On the other hand, the SHC takes a large negative 
value in Ta, Nb, W, and Mo where rid < 5. In transition metals, a conduction electron acquires the 
trajectory-dependent phase factor that originates from the atomic wavefunction. This phase factor, 
which is reminiscent of the Aharonov-Bohm phase, is the origin of the SHC in paramagnetic metals 
and that of the anomalous Hall conductivity in ferromagnetic metals. Furthermore, each transition 
metal shows huge and positive d-orbital Hall conductivity (OHC), independently of the strength of 
the spin-orbit interaction (SOI). Since the OHC is much larger than the SHC, it will be possible to 
realize a orbitronics device made of transition metals. 

PACS numbers: 72.25.Ba, 72.25.-b, 75.47.-m 



INTRODUCTION 



In recent years, spin Hall effect (SHE) in transition 
metals has received considerable attention due to its fun- 
damental as well as technological interest. The SHE is 
the phenomenon that an electric field induces a spin cur- 
rent in a transverse direction. Recent experimental ef- 
forts have revealed that many metallic compounds show 
sizable spin Hall conductivity (SHC) 0, 1 1, i . In par- 
ticular, Pt shows a huge SHC at room temperature [llH. 
The observed SHC in Pt is - 240 fie" 1 • ft" 1 cm -1 [J], 
which is about 10 4 times larger than the SHC reported 
in n-type semiconductors. This unexpected experimental 
fact cannot be understood based on the simple electron 
gas models for semiconductors 0, @. To elucidate the 
origin of the huge SHE in transition metals, several au- 
thors have studied the SHC based on the multiorbital 
tight-binding models @, HI and the first-principles band 
calculation Q. Reference [?J has shown that the orbital 
degrees of freedom in transition ions, which are absent in 
electron gas models, are crucial to realize the huge SHE 
in various transition metals. 

The SHE has a close relationship to the anomalous Hall 
effect (AHE) in the presence of the magnetization M, 
where charge current is induced by an electric field E in 
parallel to M x E. In 1954, Karplus and Luttinger (KL) 
[lfj ] solved the kinetic equation for 1"-, |-spin electrons 
in the multiband system with the z-component of the 
atomic spin-orbit interaction (SOI) Ai ^\(Z z s z );. Under 
the electric field E || y, they showed that the f-spin and 
J,-spin electrons move to the opposite direction in paral- 
lel to the x-axis. Their analysis strongly suggests that 
the AHE (SHE) occurs in ferromagnetic (paramagnetic) 



multiband systems with the SOI. The Hall effect studied 
by KL, which is due to the interband particle- hole excita- 
tion and is independent of impurity scattering, is called 
the "intrinsic Hall effect" . However, this explanation ac- 
cording to the KL-theory is too naive in that they omit- 
ted the x, y-components of the SOI A2 Y2i(lxS x + l y s y )i. 
In the case of A2 ^ 0, there is no simple relation between 
SHE and AHE since s z = ±1/2 is not a good quantum 
number. 

Note that KL did not make a mention of the SHC. The 
analogous relationship between AHE and SHE was first 
pointed out by Dyakonov and Perel [Til ]. 

After KL, theories of the intrinsic AHE [H, [lj, Q3, 

f»JSE.IISS^lJ23, HI [H, [H and the intrinsic SHE 
[fj, [j], la [24|, |25j, |2fJ have been improved based on sev- 
eral specific theoretical models. Kontani and Yamada 
[l]| studied the intrinsic AHE based on the periodic An- 
derson model by considering the SOI unperturbatively. 
Using the microscopic Fermi liquid theory, they derived 
the general expression for the anomalous Hall conduc- 
tivity (AHC) by considering all the self-energy correc- 
tion and the current vertex correction (CVC). Their 
study had clarified that the intrinsic AHE (due to the 
KL mechanism) remains finite even if all the scattering 
processes are taken into account rigorously, against the 
Smit's claim [13|. The obtained general expression has 
succeeded in explaining the huge AHC observed in heavy- 
fermion systems. It was found that the large anomalous 
velocity, which is not perpendicular to the Fermi surface 
and is the origin of the AHE, is caused by fc-derivative 
of the phase factor in c-f mixing potential. That is, the 
/-orbital degree of freedom is significant for the AHE. 
Later, AHE in (i-electron systems has been studied in- 
tensively [13, US EH, [23. 
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Recently, Murakami et al. [f| and Sinova et al. 0| 
calculated the intrinsic SHC in the Luttinger model and 
the two-dimensional Rashba model, respectively. Later, 
several authors have studied the disorder effect on SHC 
[24! l28l |29| . Inoue et al. [23[ proved that the intrinsic 
SHE in the Rashba model vanishes due to the cancella- 
tion by the CVC due to impurities. In analogy to the 
quantum charge Hall effect, it has been predicted that 
large (and quantized) SHC may be realized in massless 
Dirac electron systems, when the chemical potential lies 
inside the gap induced by the SOI. This mechanism has 
been predicted to be realized in some semiconductors [30| 
and in graphene [3l| . 

However, in usual metallic systems, existence of the 
Dirac point just at the Fermi level cannot be expected in 
general. Therefore, other novel mechanism for large SHC 
is expected to be realized in Pt and other transiton met- 
als. In fact, authors in ref. [3] presented the first report 
on the theoretical study of SHE in transition metal com- 
pound Sr2Ru04, which has no Dirac cone type dispersion 
at the Fermi level. They had found that the origin of the 
huge Hall effect is the "effective Aharonov-Bohm (AB) 
phase" induced by <i-orbital degrees of freedom, which 
are absent in semiconductors and in graphene. This 
mechanim is expected to produce huge SHE in various 
multiorbital (i-electron systems universally. Later, Ref. 
[1, Q succeeded in reproducing the SHC in Pt theoreti- 
cally. 

In the present paper, we study the intrinsic SHE in 
various 4c?- and 5c?-transition metals by taking account 
of their realistic band structures. We employ the Naval 
Research Laboratory tight-binding (NRL-TB) model [H, 
[33j , which enables us to construct nine-orbital (s + p + 
d) tight binding models for each transition metal. We 
find that both Pt (5c? 9 ) and Pd (4c? 9 ), which have face- 
centered cubic (fee) structure, show large positive SHCs. 
On the other hand, the SHCs take large negative values 
in Ta (5d 4 ), Nb (4&? 4 ), W (5d 5 ), and Mo {Ad 5 ), which 
have body-centered cubic (bec) structures. We find that 
the SHC changes smoothly with the electron number n = 
n s + rid regardless of the changes of the crystal structure, 
where n a represents the number of electrons on a-orbital. 
Among them, Pt shows the largest SHC in magnitude 
in the low resistive regime. In usual, intrinsic SHC is 
independent of resistivity in the low resistive regime (p < 
50 /iilcm), whereas it decreases in proportion to p~ 2 in 
the high resistive regime. However, we find a condition 
that the intrinsic SHC decreases as p approaches zero in 
the low resistive regime. This anomalous phenomenon 
may be realized in Ta. 

Furthermore, we study the d-orbital Hall effect (OHE), 
which is the phenomenon that an electric field induces a 
d-orbital current in a transverse direction 0, [|[ . We find 
that the d-orbital Hall conductivity (OHC) is almost one 
order of magnitude larger than the SHC, since the OHC 
occurs even in the absence of the SOI. Using the large 
OHE in transition metals, we will be able to constract a 
orbitronics device made of transition metals. In a later 



publication, we will present an intuitive (semiclassical) 
explanation for the origin of the OHC [34j . 

Finally, we comment on the extrinsic Hall effect. In 
1958, Smit [35[ studied the AHE due to the asymmetric 
scattering around the impurity in the presence of spin- 
orbit coupling, which is called the skew-scattering mech- 
anism. The AHC due to the skew-scattering is propor- 
tional to p^ 1 if elastic scattering is dominant. In 1970, 
Berger proposed another mechanism of extrinsic Hall ef- 
fect, the side jump due to impurities [36|. This mecha- 
nism gives the AHC in proportion to p~ 2 . Both extrinsic 
Hall effects vanish where the inelastic scattering due to 
electron-electron or electron-phonon interaction is dom- 
inant over the elastic scattering. Both mechanisms (the 
skew-scattering and the side jump) cause the extrinsic 
SHE [U [H, [H, H^. In the present paper, we do not 
study the extrinsic SHE, which is sensitive to the char- 
acter of the impurity potential. It is an important future 
problem to study the extrinsic SHE in realistic multior- 
bital tight-binding models. 



II. MODEL AND HAMILTONIAN 

In the present study, we use the NRL-TB model 
[32I [33j to obtain the band structure in various transition 
metals. Here, we shortly explain this model. The NRL- 
TB model employs the scheme of two-center and non- 
orthogonal Slater-Koster (SK) Hamiltonian 0. The SK 
parameters are represented with distance- and environ- 
ment dependent parameters that are determined so that 
the total energy and the band structures agree with those 
obtained by full-potential LAPW LDA calculations. The 
root mean square error in the fitting is about 0.002-0.004 
Ry [12] • This fitting error is small enough to perform re- 
liable numerical calculations. To describe the electronic 
state in Ad (5c?) metals, we consider 5s, bp, and Ad (6s, 
6p, and 5d) orbitals, that is, we consider nine orbitals 
per atom. Hopping integrals between a pair of atoms are 
then expressed with ten SK parameters (ssa, spa, ppu, 
pprr, sda, pda , pdn, dda, ddn, and ddS). In this study, we 
consider hopping (and overlap) integrals between up to 
sixth nearest neighbor sites for metals with fee and bec 
structures, and ninth nearest neighbor sites for hexagonal 
closed packed (hep) structures. 

Table [J shows the crystal structure, electron number 
per atom, and the coupling constant A of SOI A Yl - Ijij 
(i = x,y,z) for various Ad- and 5c?- transition metals. 
(jridd nd m s s ns ) represent the electronic configuration of 
an isolated atom, where rrid and m s is the main quantum 
number, and n s and rid is the number of electrons on s- 
and G?-orbital. In this table, we put A=0.03 Ry for 5c? elec- 
tron in Pt and A=0.013 Ry for 4g? electron in Pd, accord- 
ing to optical spectroscopy [43| . For other 4c?- and 5d- 
transition metals, we used Herman-Skillman atomic spin- 
orbit parameters (44|: These parameters had been cal- 
culated by using the self-consistent Hartree-Fock-Slater 
atomic functions. Here, we consider only the d-orbital 
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SOI, and neglect other SOI terms which may possess a 
fc-dependence. Hereafter, we set the unit of energy Ry; 
1 Ry = 13.6 eV. 



TABLE I: The crystal structure, electron number per atom, 
and the coupling constant A of SOI for various transition met- 
als. (mdd nd rn s s nB ) represent the electronic configration of an 
isolated atom. Here, md and m s is the main quantum number, 
and n a and rid is the number of electrons on s- and d-orbital. 
Bcc, hep and fee represents a body-centered cubic, hexagonal 
closed packed and face-centered cubic, respectively. 



metals 



structure 



electron number 



SOI (Ry) 



Nb 


bcc 


5 (4d 4 5s 1 ) 


0.006 


Mo 


bcc 


6 {4d B 5s 1 ) 


0.007 


Tc 


hep 


7 (4d 6 5s 1 ) 

8 (Ad^s 1 ) 


0.009 


Ru 


hep 


0.01 


Rh 


fee 


9 (Ad 8 5s 1 ) 


0.011 


Pel 


fee 


10 (4d 10 5s°) 


0.013 


Ag 


fee 


11 (4d 10 5s 1 ) 


0.019 


Ta 


bcc 


5 (5d a &s z ) 

6 (5d 4 6s 2 ) 

7 (5d 5 6s 2 ) 


0.023 


W 


bcc 


0.027 


Re 


hep 


0.025 


Os 


hep 


8 (5d 6 6s 2 ) 


0.025 


Ir 


fee 


9 (5d 9 6s°) 


0.025 


Pt 


fee 


10 ^d^s 1 ) 


0.03 


Au 


fee 


11 (5d 10 6s 1 ) 


0.03 



In the presence of SOI for Ad or 5c? electrons, the total 
Hamiltonian is given by 



H 



( H + Xl z /2 
V X(L + il y )/2 



X(l x - ily)/2 

Ha — Xl z /2 



(1) 



where the first and the second rows (columns) correspond 
to s z = +h/2 (f-spin) and s z = —h/2 Q-spin). Hq is a 
9x9 matrix given by NRL-TB model for bcc and fee 
structure. In the case of hep structure, Hq is an 18x18 
matrix, since a unit cell contains two atoms. The matrix 
elements of I for d-orbital are given by [43| 
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where the first to the fifth rows (columns) correspond to 
e?-orbitals xy, yz, zx, x 2 — y 2 , and 3z 2 — r 2 , respectively. 

In NRL-TB model, we use the non-orthogonal basis 
since the atomic wave functions of different sites are not 
orthogonal: 



dr4>* a {r - Ri)<t>p{r - R«) = O aP (R t - Ri>), (5) 



where <j) a {r — Ri) represents the atomic wave function 
at the itli site and a, f3 is orbital state indices, and 
O a p{Ri — Ri') represents the overlap integral between 
different sites. When the overlap integrals between dif- 
ferent sites are negligible, eq. ([5]) is simplified as 



O a p{Ri — Ri') — 5 a p5ni 



(6) 



This approximation is rather appropriate when a and 
(3 correspond to d-orbitals, since d-orbital atomic wave 
functions are localized well. However, when a and [3 are 
either s- or p-orbital, O a p(Ri — Ri/) is large even when 
i 7^ i', and therefore eq. ([6]) is not satisfied. 

The band structures obtained for the present model in 
Au and Ta are shown in Fig. [1] These band structures 
are derived by taking the overlap integrals between dif- 
ferent sites into account correctly. The methods how to 
calculate these band structures are explained in section 
MI CI Near the Fermi level, the obtained band structures 
are in good agreement with the results of the relativistic 
first- principles calculations [H, |46[ . Since the Fermi sur- 
face is mainly composed of d-electrons in the transition 
metals, in which we study the SHE and OHE, the band 
structure near the Fermi level is described well when 
O a f3{Ri — Rf) is approximated by eq. ([6]). Therefore, 
the calculations of SHE and OHE using eq. © seems to 
give semiquantitatively reliable results [8j. However, for 
a more quantitative study of SHE and OHE, we need to 
consider the overlap integrals between different sites. 

Until section MI C| we will use the simplified overlap 
integrals given by eq. © to simplify the explanation. In 
section ITlI CI we will study the SHE and OHE by con- 
sidering the overlap integrals in eq. ([5]) correctly. The 
18x18 matrix form of the Green function without impu- 
rities is given by G°(k, u) = (uj + [i — H)~ x where \x rep- 
resents the chemical potential. There is a fc-dependent 
unitary matrix U which diagonalize the Hamiltonian H 
as follows: 



a/3 



U} a H a Pu 0m = E l k S, 



lm •) 



(7) 



(3) 



where a, (3 is the orbital indices, and I and m are the 
band indices. 
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FIG. 1: (upper panel) Band structure of Au. Here, 
r = (0,0,0), X=(tt,0,0), W=(7r l7 r/2,0), L=(7r/2,ir/2,7r/2), 
and K=(37r/4, 37r/4, 0). (lower panel) Band structure of 
Ta. Here, F = (0,0,0), H=(tt,0,0), N=(tt/2, tt/2, 0), 
P=(7r/2, 7r/2, 7r/2). Near the Fermi level, we see that the 
band structures obtained in the present model agree well with 
the result of the relativistic first-principles calculation in refs. 
11. 



Here, we consider the quasiparticle damping rate 



r, which is given by 



Sfc(-iO) - £ fc (+«0) /2i, where 



T.k(uj) is the self-energy matrix. In the Born ap- 
proximation for / <C Wband, f is given by V = 
nim P I 2 ^i [g(—i0) — g(+i0)], where n imp is the impurity 
concentration, / is the impurity potential, Wband is 
the bandwidth, and g(w) is the local Green function 
g(u) — j^^2 k G{k 1 uo) 1 respectively. Here, N is the to- 
tal number of lattice points. In the T-matrix approx- 
imation for general strengh of /, T is given by T — 

n imT) I — { 1. The retarded 

mp 2i\l-Ig(-i0) l-lg(+i0)j 

and advanced Green functions are given by 

G R (k,uj) = (uj +fi-H + it)- 1 , 
G A (k,oj) ^(uj + fi-H-iT)- 1 . 



(8) 
(9) 



In the present model, the charge current operator for 
/i-direction (/i = x,y) is given by 



1° 







o £ 



(10) 



where — — eff 2 -, and — e (e > 0) is the elec- 
tron charge. Here, atomic SOI is not involved in 
the charge current operator since it is fe-independent. 

Also, the s 2 -spin current operator — \ J^,s 



{^ C A}/2 



(jc§ z + s z jf\ /2 [13,|43 is gi ven b Y 



J S u = (-n/e) 



o 



(11) 



o -# 

and Z z -orbital current operator is given by 

4° = {4 c X}/2. (12) 

We will discuss the validity of these current operators in 
Appendix [B] in more detail. 



III. SPIN AND ORBITAL HALL 
CONDUCTIVITIES 



SHC and OHC without overlap integrals 
between different sites 



In this section, we derive the general expressions for 
the intrinsic SHE and OHE based on the linear-response 
theory. As we will discuss in IIII B[ we can safely neglect 
the CVC in calculating SHC and OHC in the present 
model. Therefore, the SHC at T=0 is given by a z xy — 
according to Streda [48[ , where 



zl , zll 
xy 1 xy 



r zII 



1 



2ttjV ^ 

k 



Tr \J°G 



S/Vi? jC 



jy G 



j=0 



(13) 



AttN 



E 



cLuTt 



- J„ G J„ 



C dG R 



A) 



(14) 



Here, / and II represent the "Fermi surface term" and 
the "Fermi sea term" , respectively. In the same way, the 
OHC of the Fermi surface term Oil, and that of the Fermi 
sea term 0*y are respectively given by eqs. (fT3| and (fl4|) 
by replacing J x with the / z -orbital current operator J® 
in eq. p^|). 

In the Born approximation, the quasiparticle damping 
rate depends on orbital indices. If A = 0, T is diagonal 
with respect to orbital: T a p = 7 Q <5 a ,3, where a and (3 are 
orbital indices. Then, the offdiagonal terms are negligible 
when A <C Wband- la is the quasiparticle damping rate 
for a-orbital, and is proportional to the local density of 
states (LDOS) for a-orbital, p a (Q). In the present model, 
d-orbital LDOSs are about equal in magnitude. There- 
fore, 7 Q is approximately independent of a and can be 
approximated by a constant 7: T a p = 7<5 Q , j a. Using this 
constant 7 approximation, we derive the general expres- 
sions for the Fermi surface term and the Fermi sea term. 
In this case, the retarded and advanced Green function 
can be diagonalized using unitary matrix U given by eq. 
([7]) as follows: 



J2 U L°: 



It 



Ufim — 



>lm 



a/3 



(15) 
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where I, mis the band indicies. Therefore, we can rewrite 
eqs. (|13l) and ([14")) by using U as follows: 



z/ = 



^ E (■tfrWV 



i 



2ttN 



k,l^rn 



2irN 



1- £ Im{(^r(J^} 



xlm 



(^-»7)( J E2 l +i 7 ) 



(16) 



27riV 



E 

k,l^m 



xlm 



dcjlm{(J^ ml (J^) lm } 

1 



1 



(w - + 17) (w - + 27 ) 2 



(17) 



where (Jf ) ml is given by £ Q/3 U^J^U^. Note that 
we dropped the diagonal terms I = m in the summa- 
tions in eqs. (fTr?)) and (fTTj) since they vanish identi- 
cally. We also note that the transformation from the 
first row to the second row in eq. (| 16(1 was performed 
since J2i m R e {(Jx) ml (Jy) lm } vanishes identically after 
fc-summation. After performing the cj-mtegration in eq. 
(fl7|) . the Fermi sea term is given by a^ 11 = cr** Ia + a zI y Ib , 
where 



a 



■<:y 



2ttN 



E l k + EZ-2i~t 



X Im 



(4-i 7 )(^-i 7 ) 



(18) 



zllb 



E l k -il 



k,l^m 

x Im On 



(4 - e%) 2 



(19) 



Here, we used the following relation to perform the u- 
integration: 

dx 



(20) 



1 



1 



(x — a) 2 (x — b) (x — a)(x — b) 2 
a + b — 2/i 2 / a — p, 



rln 



(a — b)(a — ji){b — p) (a — b) 2 \b — p, 



In the case of 7 — > 0, cr^f/ 6 given by eq. corre- 
sponds to the Berry cuvature term given by: 



< /b = ^E/(^( fc )' 



(21) 



where tt l (k) represents the Berry curvature given by 



o*(fc) = E 



2Im 



(K - E™) 2 



(22) 



The relation a* = a z i Ib has been frequently assumed 
in literatures such as ref. [9J. In this study, we calcu- 
late all the Fermi surface and Fermi sea term correctly, 
and elucidate how each term contribute to the SHC and 
OHC. In section IVBl we will discuss the 7-dependece of 
v z Jy,a z J y la and a zI y Ib in detail. 



B. Discussion on the CVC 



In section UlI A[ we have neglected the CVC. Here, we 
calculate the CVC due to the local impurity potential in 
the Born approximation, and show that it is negligible in 
transition metals. In the Born approximation, the lowest 
order CVC is given by 



1 

N 



'■imp 



i 2 J2g r J^'g' 



(23) 



Its diagrammatic expression is given in Fig. O The mag- 
nitude of CVC depends on the model. For example, the 
CVC vanishes identically in the d-orbital tight-binding 
models with atomic SOI @, H3l- I n contrast, the CVC 
plays an essensial role in a Rashba model: the SHC van- 
ishes due to the cancellation by CVC [23j. 

Here, we study the CVC in fee and bec transition 
metals, where each atomic site is a center of inversion 
symmetry. The s- and d- orbital atomic wave functions 
have an even parity with respect to k — > — k, whereas 
the p-orbital atomic wave functions have an odd parity. 
Therefore, the Hamiltonian for a fee and bec metal has 
a following relationship: 



(H (k)) a p = p(H (-k)) a p, 



(24) 



where p = — 1 only when either a or (3 is p-orbital, oth- 
erwise p=l. It is easy to show that (G(k, uo)) a i3 has the 
same parity with (Ho(k)) a p. Therefore, when both a 
and /3 are (s,d)- orbitals, (AJ^) a p = since the {a, (3)- 
components of (d/dk^G = GJ^G is an odd function. 
On the other hand, (A J^) a p 7^ when either a or (3 
is p-orbital since (dG / dk^) a p is an even function. Al- 
though (AjC) a p originating from the p-orbital is finite, 
it is small in magnitude since the bp, 6p-level is about 
20eV higher than the Fermi level p and the p-electron 
density of states (DOS) at p, is very small in all transi- 
tion metals. 



FIG. 2: The diagrammatic expressions of the current vertex 
correction due to the local impurity potentials. Here, the 
diagram represents the lowest order correction to the current 
Aj£ and spin Hall conductivity Aa% y . 
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Now, let us verify numerically that the contribution of 
CVC to SHC is small in magnitude. The correction to 
the SHC due to the lowest order CVC is given by 

k,k' 

- 1 r> t 2 \^Tt \r A f s n R n R i c n A 

— 2ttN 2 mp Z-' i y ' 

(25) 

Its diagrammatic expression is given in Fig. [2] We cal- 
culated Act^j, numerically and found that it is very small 
compared to o z xy without CVC: the ratio | Aa* y \ / 
| a z xy | is - 0.02 for Ta and Pt, and - 0.005 for W 
when 7=0.002^0.02. The ratio is independent of 7 by 
the following reason: Since ^ k G R G A ~ 0(7 _1 ), and 
7 oc n,i m pl 2 , Ajff in eq. ([2"3"| is independent of 7. In the 
present model, the higher order correction to the SHC 
should be negligible. 

In the case of hep transition metals, (AJf) a/ g ^ even 
when a and (3 are (s, <i)-orbitals, since each atomic site 
is not a center of inversion symmetry. To find out the 
importance of the CVC in hep metals, we performed the 
numerical calculation for Os, and found that the ratio 
I I / I a% y I is ~ 0.06 in Os in the low resistive 

regime. Although it is much larger than that in Pt, Ta 
and W, the CVC is qualitatively negligible even in hep 
transition metals. Therefore, we are allowed to neglect 
the CVC even for hep transition metals. 

C. SHC and OHC considering the overlap integrals 
between different sites 

In the previous section, we studied the SHC and OHC 
under the assumption that the atomic wave functions of 
different sites are orthogonal. However, for a more accu- 
rate quantitative study of the intrinsic SHE and OHE, 
we need to take the off-diagonal elements of O a p in eq. 
(J5]) into account correctly. In this subsection, we explain 
how to calculate the SHE and OHE when the overlap 
integrals between different sites are considered. 

Here, we introduce the Bloch wave function which is 
given by 

<f> ka (r) = -j=Y^e ih - R *<f> a (r-R i ). (26) 

v i 

In this case, the inner product between the Bloch wave 
functions with different k and a is given by 

J drcf>la(r)4>h^(r) = 6 kk ,O a(3 (k), (27) 

where 

o aP (k) = Y. c ~ tkRi °^ R *)- ( 28 ) 



Here, O a p(Ri) in eq. I)28p is the overlap integral which is 
defined by eq. ([5]). Therefore, when the overlap integrals 
between different sites are considered, the Bloch wave 
function given by eq. (|26p is non-orthogonal. 

By including the chemical potential /x, the kinetic term 
of Hamiltonian is given by [491 ] 

Ho=J2 C L [hap(k) - fxO a p(k)] c k0 , (29) 

h,a,P 

where c ka is defined by 

c fcQ = -L5> lfcR *c lct . (30) 

v i 

Here, Cj a is an annihilation opertaor of an electron in the 
a orbital state at ith site. As the atomic wave functions 
at different sites are non-orthogonal, creation and anni- 
hilation operators ct , c ka do not satisfy the canonical 
anticommutation relations, but instead satify (49| 

{cka,cl, } = Skk>0~p(k). (31) 

Since matrix O(fc) is positive definite Hermitian matrix, 
we can introduce a matrix S(k) that transforms 0(k) 
into the unit matrix 1: 

&(k)6(k)S{k) = 1. (32) 

We note that matrix S(k) cannot be determined 
uniquely: Using an arbitrary unitary matrix X, S' = SX 
also satisfies eq. (|32[) . Here, we introduce the following 
new basis (cka,c ka ) using S(k): 

Cka = ^^2s~^(k)c kf j. (33) 

P 

We can easily verify that these operators (c ka ,c ka ) 
satisfy the canonical anticommutation relations 
{cfea,Cfe //3 } = 6 kk '5 a p. In this basis, eq. is 
rewritten as 

Hp = ^ 4a \hai3{k) ~ M<W] Cfc/3, (34) 

where h a p{k) = {S\k)h{k)S{k)) af3 = (h(k)) a ( } . There- 
fore, the Green function in the (c ka ,c k ) basis is given 
by 

G(k, u>)=(w + (i- h{k)^j 1 . (35) 

Next, we derive the expression for the current opera- 
tor in the (c ka , c ka ) basis. From the continuity equation 
■§^n(r) + V • j(r) — 0, we obtain 

d 

—n(q) = -iq-j(q), (36) 



7 



where n is electron number density, -^n(q) can be calcu- 
lated by the equation of motion: -^n(q) — i[H,n(q)]. 
Therefore, x-component of the current operator j x is 
given by 



(37) 



10 1 



j x = Urn [H, n(q) 



In Appendix 13 we show that n(q) is given by 

n{q) = °^( k ) c {- q /2, a c k+ q /2j3, (38) 

k,a,/3 

which is an exact expression for the first order of |q|. Us- 
ing the following relationship: [AB, CD] = A{B, C}D - 
AC{B, D} + {A, C}DB - C{A, D}B, [H, n(q)} is given 
by 

[H,n(q)} = A- q /2, a mk-q/2)6-\k-q/2)d{k)rP 

k,a,0 



without overlap 

o_o £ „ lOOOfee-^-icm- 1 
o-o- : -o-o - . , 

■ > .7 

— sl Q -. 




10" 2 
10 



10- 2 10" 1 
y (Ry) 



{6{k)6- l {k + q/2)h(k + q/2)) Q/3 }c fc 



+<?/2,/3- 

(39) 



By substituting eq. (|3"9")) into eq. (|3"Tj) . we obtain the ex- 
pression for the velocity in the (cfe Q , ct ) basis as follows: 

ix[k) = ^fc7 + + \DUk)h(k), (40) 

where D x {k) is given by 

Mk) = {^-d- 1 (k)}d(k) = -d- 1 (k)^-d(k). 



(41) 



Apparently, D x (k)—Q in an orthogonal basis. We call the 
second and the third terms in eq. (|40|) the overlap inte- 
gral current. In the (cfc a , c k ) basis, the velocity v x (k) is 
given by 



v x (k) = &(k)% x {k)S{k). 



(42) 



Therefore, even when the overlap integrals between dif- 
ferent sites exist, we can calculate the SHC and OHC in 
the basis (5ka,c ko ,) using the matrix S(k). In this basis, 
Green function G(k,oj) is given by eq. ([33)1 . The charge 
current operator J~j and the spin current operator 
are given by 



FIG. 3: 7-dependence of SHC in Pt obtained by using eq. 
© (with overlap) and that obtained by using eq. (0 (with- 
out overlap). When the overlap integrals between different 
sites are considered, the magnitude of SHC is reduced by 
about half. The resistivity that corresponds to 7 = 0.01 is 
~ 160/if2cm in Pt. 



Figure [3] shows the obtained SHCs in Pt by considering 
the overlap integrals between different sites in NRL-TB 
model. For comparison, SHC given by using eq. ([5]) is 
also shown. We find that the magnitude of the SHC is 
reduced by about half when the overlap integral is con- 
sidered correctly. We have verified that this is mainly due 
to the changes of the band spectra, whereas the modifica- 
tion due to the overlap integral current, which is given by 
the second and third terms in eq. ptj]). is less than 10% 
in magnitude. Since the CVC is also little affected by the 
modification of the velocity, we can also safely neglect the 
CVC when the overlap integrals between different sites 
are considered. 

Here, we comment on the previous study ref. 0]. 
Therein, they studied the SHC in Pt using eq. Q . Since 
the Fermi surface is mainly composed of d-electrons in 
Pt and the <i-orbital atomic wave functions are well local- 
ized, the band structure near the Fermi level is described 
well in this approximation. Therefore, the calculations 
of SHC and OHC in the absence of the overlap integrals 
between different sites give semiquantitatively reliable re- 
sults, which can be recognized from Fig. [3J 



IV. NUMERICAL STUDY 



A. Numerical results for SHC and OHC 



J C ,= 



Vu 



l 7< 







-v,j 
"o _ 



(43) 



Also, the ^-orbital current is given by J® = {J^,l z }/2 

where l z = &(k%S(k). Therefore, SHC and OHC can 
be calculated by substituting G(k, to), j£ and in eqs. 

(IT31) and (HP with G(k,w), J° and Jj*, respectively. 



In this section, we perform the numerical calculations 
for the SHC and OHC in various Ad- and 5d- transition 
metals by considering the overlap integrals between dif- 
ferent sites given by eq. ([5]). In particular, we clarify SOI 
(A), the quasiparticle damping rate (7) and the chemical 
potential (/x) dependences of the SHC and OHC in each 
metal. Here, we note that the unit of the SHC and OHC 
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FIG. 4: A-dependence of SHC and OHC in Pt given by the 
Born approximation and constant 7 approximation. The 
SHCs obtained in these two approximations give quantita- 
tively similar results. We stress that the OHC is finite even 
if A = 0. 

First, we discuss the quasiparticle damping rate T de- 
pendence of SHC. In the Born approximation, T depends 
on orbital index. When A <C Wb an d, T is diagonal with 
respect to orbital: T a p — 7 Q <5 Q /3, where j a oc p a (0). 
On the other hand, the quasiparticle damping rate j a 
is indepedent of orbital in the constant 7 approximation: 
T a/3 = jS al3 . In Fig. |U the SHCs for 7 = 0.007 in 
these two approximations are shown. We see that the 
obtained SHC is quantitatively similar in both approxi- 
mations. This fact can be explained as follows: In tran- 
sition metals, the LDOS of ti g - (d xy ,d yz ,d zx ) and e g - 
(d x 2 _ y 2 , d 3z 2 _ r 2) orbitals are about equal in magnitude. 
Since 7 Q oc p a (0) in the Born approximation, two ap- 
proximation give the similar results in transition metals. 
For this reason, we use the constant 7 approximation 
hereafter. In contrast, the SHC in Sr2Ru04 given by the 
Born approximation is much larger than that given by 
the constant 7 approximation since a-dependence of p a 
is large Q. 

Figure 2] also shows the OHCs in these two approxi- 
mations. As already pointed out in Refs. @, 3, a huge 
OHC appears even if A = 0, and it slowly increases with 
A. In Sr 2 Ru04, in contrast, the OHC slowly decreases 
with A 0- 

Figure [5] shows the electron number n-dependence of 
the SHC, where n = n s + n^. Note that the crystal 
structure of various transition metals is shown in Table 
HI The SHC obtained in the present model for 7 = 0.002 
is shown in Fig. [5] (a), and for 7 = 0.02 and 0.2 in Fig. 
[5] (b). The SHC is negative for n = 5, 6, and positive for 
n = 9 ~ 11: The SHC changes its sign around n = 7 
and 8. The magnitude of SHC is largest in Pt for 7 = 
0.002 and 0.02, where the corresponding resistivities are 
~ 8pilcm and ~ 64/zf2cm in Pt, respectively. When 7 = 
0.2, however, the absolute values of SHCs in Ta and W 
become larger than that of Pt, where p ~ 220 (250) pilcm 




W -0.0$ 4d 
-0. 



J I I 1_ 



n 



10 11 



4d Nb Mo Tc Ru Rh Pd Ag 
5d Ta W Re Os Ir Pt Au 



FIG. 5: n-dependence of SHC for 7 = 0.002, 0.02 and 0.2. In 
(a), we see that Pt shows the largest SHC in magnitude for 
7 = 0.002. The open symbols represents the SHC in Pt for 
n — 5 ~ 9. In (b), the SHCs obtained in the present model 
for n — 7,8 (hep structure) are also shown. SHC in W takes 
the largest value for 7 = 0.2. 



in Pt (Ta and W). Therefore, large negative values of 
SHCs in Ta and W will be observed even in high resistive 
samples. For comparison, we also calculated the SHC 
for n = 5 ~ 9 using the band structure of Pt, which is 
represented as the open symbols in Fig. [5] (a) . We see 
that the magnitude of SHC in Pt with n = 9 does not 
reproduce that in Ir. The same is true in Ta and W (bec 
structure). Therefore, we need to calculate SHC using a 
correct band structure for each metal. 

Next, we examine the A-dependence of the SHC. We 
verified that the SHC in each metal increases approx- 
imately proportional to A as shown in Fig. [5] To 
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elucidate the origin of SHC, we calculated the SHC 
when SOI is anisotropic: Ai ^2{l z s z ) + A2 ^2{l x s x + lySy)- 
We find that the SHC for H so = A£(f z s z ) (Ai = 
A, A2 = 0) is as large as the SHC in the isotropic case 
(Ai = A2 = A). In contrast, the SHC for Hso = 
\^2{l x s x + lySy) (Ai = 0, A 2 = A) is one order of mag- 
nitude smaller than the isotropic case. Therefore, the 
^-component of the SOI gives the decisive contribution 
to the SHC. The matrix element of l z is finite only for 
(yz\l z \zx) = -(zx\l z \yz) = i and (xy\l z \x 2 - y 2 ) = 



(x 2 - y 2 \l z \xy) = 2i. Note that d 



XI) 



and d x 2_ y 2 



orbitals (d yz - and dzz-orbitals) are given by the linear 
combination of l z = ±2 (l z = ±1). Here, we examine 
which orbitals among them cause the significant contri- 
bution to the SHC. z-component of SOI is rewritten as 

A 3 E { P{ £ = l ) fa) } t + A 4 E = 4 ) fa') 

i i 

where P(l 2 — n) represents the projection operator. SHC 
caused by d xy - and c? x 2_ J( 2-orbitals is given by setting 
A3 = 0, A4 = A, which is represented as l z = ±2 in 
Table HU Similarly, SHC caused d yz - and d za; -orbitals is 
given by setting A3 = A, A4 = 0, which is represented 
as l z = ±1 in Table [TTJ We see that the interorbital 
transition between d xy - and d x i_ y i- orbitals causes the 
significant contribution to the SHC in many metals. Only 
in the case of Mo, W and Ir, the contribution of d zx - and 
oL r -orbitals is comparable to that of cL„- and d, 



-y 



orbitals. In other metals, d xy - and d x i_ y i- orbitals give 
the dominant contribution to the SHC. 



d x 2 _ y i -orbitals. Here, we set 7 = 0.02. l z = ±2 (l z = ±1) 
represents the SHC caused by d xy - and d x 2_ y 2- orbitals (d yz - 
and d zx - orbitals). The ratio represents (SHC from l z — ±2) 
/ (SHC from l z = ±1). We see that d xy - and d x 2_ y 2- orbitals 
cause the significant contribution to the SHC in many metals. 



metals 


l x = ±1 


l z = ±2 


ratio 


Nb(4d 4 5s x ) 
Mo(4d 5 5s 1 ) 
Rh(4d 8 5s 1 ) 
Pd( 4d 10 5s°) 
Ag(4d 10 5s 1 ) 


-0.0332 
-0.0474 
0.0847 
0.0847 
0.00224 


-0.0770 
-0.0587 
0.269 
0.455 
0.0181 


2.32 
1.24 
3.18 
5.37 
8.08 


W(5d 4 6s 2 ) 
Ir(5d 9 6s°) 
Pt(5d 9 6s 1 ) 
Au(5d 10 6s 1 ) 


-0.0222 
-0.174 
0.0123 
0.136 
0.0177 


-0.254 
-0.205 
0.0231 
0.678 
0.0987 


11.4 
1.18 
1.89 
4.98 
5.59 



Here, we show the OHCs for 7 = 0.02 in various tran- 
sition metals in Fig. [6] We see that all the Ad and 5d 
transition metals show huge and positive OHCs, which 
is almost one order of magnitude larger than that of the 
SHCs. In Au (Ag), the OHC takes a small value since 
the d-electron DOS is small at the Fermi level. Therefore, 
huge and positive OHC is a universal nature of transition 
metals. As in the case of the SHE, the intrinsic OHE 
shows the crossover behavior: the OHC independent of 




FIG. 6: n-dependence of OHC for 7 = 0.02. The obtained 
OHCs are positive for all metals, and they are about 10 times 
larger than the SHCs except for Pt and Pd. 



p in the low resistive regime, whereas it decreases in pro- 
portion to p~ 2 in the high resistive regime @, H|. In a 
later publication, we will present an intuitive (semiclas- 
sical) explanation for the origin of the OHC [34| . 

Now, we discuss the 7-dependences of SHC and OHC. 
The 7-dependence of intrinsic SHCs in Ta and W are 
shown in Fig. [71 In usual, the intrinsic Hall conduc- 
tivities are independent of 7 in the low resistive regime 
where 7 -C A, whereas it decreases approximately in pro- 
portion to p~ 2 in the high resistive regime where A <C 7 
[151 122I]: Here, A represents the band-splitting measured 
from the Fermi level. In W, A is ~ 0.04, and in Pt, A 
is ~ 0.035 @. We find that W shows a typical coherent- 
incoherent crossover at 7 ~ A in Fig. [3(b). We have also 
verified that the coherent-incoherent crossover behavior 
is universally seen in many transition metals including 
Pt, which is shown in Fig. [3] (a). However, as shown in 
Fig. [7j (a), SHC in Ta shows an exceptional behavior: It 
takes a maximum value around 7 ~ 0.02 and decreases 
as 7 decreases in the low resistive regime. We find that 
this anomalous behavior can arise when almost degener- 
ate anticrossing points exist slightly away from the Fermi 
level. We will discuss the reason in detail in section llVBI 
In fig. [7j (b), the 7-dependences of the Fermi sur- 
face term a z J y and the Fermi sea terms <7 zI y Ia ,a zI y Ih are 
also shown. In the low resistive regime, the relation 
a* y — & xy Ib holds well^and o z J y Ih reproduces the total 
Hall conductivity a z xy 
however, u%„ ~ a z y 



In the high resistive regime, 
, y is quite different from 
cr X y Ib in the high resistive regime. As a result, the rela- 
tionship 



■nj 



a z ~ a zI (Fermi surface term) , 



(44) 



is recognized for a wide range of 7. We will discuss the 
crossover behavior of the intrinsic Hall conductivity in 
more detail in section IV Bl 
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FIG. 7: 7-depnedence of SHC and OHC in Ta (a) and W (b). 
The corresponding resistivity p to 7 = 0.1 is ~ 190 /iflcm in 
Ta and ~ 220 /ificm in W. 



We also discuss the 7-dependence of OHC. The ob- 
tained OHCs in Ta and W are shown in Fig. [71 The 
coherent-incoherent crossover behavior of OHC is recog- 
nized in Fig. [7] In contrast to the 7-dependence of SHC 
in Ta, OHC shows a typical crossover behavior. We also 
verified that OHC is finite even if A = since the (i-orbital 
current in eq. (|12p is independent of the spin index Q. 

Now, we discuss the //-dependence of SHC in Pt and 
Ta assuming that the bandstructure is rigid. Experi- 
mentally, the chemical potential /t can be controlled by 
composing alloys. Figure [5] (a) shows the /^-dependence 
of af y Ib in Pt for 7 = 0.002,0.007 and 0.02. Here, the 
chemical potential is given by /i= /io+ A/i, where /io rep- 
resents the true value of chemical potential. We see that 
the SHC shows a peak around A/t = and it decreases 
when /i is raised or lowered from its true value. This 
A/x-dependence of SHC obtained in the present model in 
Pt seems to be in good agreement with that in ref. @. 
We see that SHC for 7 = 0.002 is about 45% larger than 
that for 7=0.02 at A/t=0. 

Here, we elucidate from which part of the surface 
the SHC in Pt originates by calculating o z xy {k) 




0.0 0.05 
An.(Ry) 




0.0 0.05 
A[i (Ry) 

FIG. 8: /i-dependence of af y Ib in (a) Pt and (b) Ta for 
7 = 0.002,0.007,0.02. The chemical potential is given by 
/i=/to+A/i where /io represents the true value of chemical po- 
tential. Total electron number n is shown in the upper hori- 
zontal axis. Note that the sign of the SHC in Ta is negative. 




FIG. 9: The ^-dependence of E l k for (k y , k z ) = (0.20?r, 0.21tt) 
in Pt is shown. The bandsplitting measured from the 
Fermi level is A ~ 0.035. A wide area around 
(k x , k y , fc z )=(0.547r, 0.20-7T, 0.2l7r) gives a dominant contribu- 
tion to SHC in Pt. 



1 E 

8 ,<V , 

k x ,k v ,k z 



Tr 



Note that 



k'w=0 



only on the Fermi surface, and it takes a huge value 
at (0.76?r,0,0) (on T-X) and at (0.447r,0.447r,0.447r) (on 
L-r) since two bands are very close near the Fermi 
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(fe) 



3000 for 



level in the present model [8j 

the former point and <J xy (k) ~ 5000 for the latter 
point in the present model. However, the contribu- 
tion of these points to the SHC is small after taking k- 
summation using 512 3 fc-meshes. We verified that the 
dominant contribution comes from a wide area around 
(k x , k y , ^)=(0.54?r, 0.20tt, 0.21tt) as shown in Fig. M 

Figure [5] (b) shows the /^-dependence of cr xy Ib in Ta for 
7 = 0.002, 0.007 and 0.02. In Ta, several "steady points", 
where the magnitude of SHC is approximately indepen- 
dent of 7 in the low resistive regime, are recognized at 
Afi = -0.02 and A/x = +0.05 in Fig. M 0). When /x 
decreases across the steady point at A/i = —0.02, the 
magnitude of SHC increases and reaches a peak around 
A/i = —0.04. This peak originates from the almost de- 
generate anticrossing bands, which is discussed in more 
detail in section HVB1 

Finally, we explain why SHC in Ir is small in mag- 
nitude, by analysing ^2 k o xy (k). Here, we divide the 
fc-summation into fc + -region and fc _ -region, where fc + - 
region (fc~-region) represents the region where <J xy (k) > 
{(J z xv {k) < 0) holds: 



k+ 



(fe) 



(fe) 



(45) 



In many transition metals, such as Pt and Ta, either cr xy ~ 



is much larger than the other. In Ir, however, 

holds, and 



or , 

we have verified that the relation a xy ~ \a xy 
therefore <rf„, becomes small in magnitude. 

Ay <-> 



B. Mechanism of impurity assisted SHC 
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FIG. 10: (a) Anticrossing bands in NRL-TB model 
for Ta. We find an accidental degenerate point at 
(kx, k y , k z )=(0, 0.12-7T, 0.337t) in the present model with A=0. 
The bandsplitting at fe* is ~ 0.003. (b) Band structure around 
an accidental degenerate point slightly away from the Fermi 
level. The region where fc-summation is performed for 7^0 
and finite 7 are represented by using (A), (B) and (C). A 
represents the band-splitting measured from the Fermi level 
and fc* represents the point of the minimum band-splitting 
around the accidental degenerate point. 



In the previous section, we have verified that the SHC 
in all Ad- and bd- transition metals except for Ta are 
independent of 7 in the low resistive regime (7 <C A). 
Here, we show that the SHC can show nonmonotonic 7 
dependence in the low resistive regime when almost de- 
generate anticrossing points exist slightly away from the 
Fermi level, due to the impurity-assisted intcrband exci- 
tation. This is the origin of the anomalous 7-dependence 
of the SHC in Ta for 7 < 0.02 in Fig. [7] (a). We call this 
phenomenon the impurity- assisted SHE. 

In Ta, there are several accidental degenerate points 
with A=0 slightly away from the Fermi level. We show 
the anticrossing bands of NRL-TB model for Ta in 
Fig. [10] (a). We find an accidental degenerate point at 
(k x , k v , k z )=(0, 0.127T, 0.337r) in the present model with 
A=0. Note that this degeneracy splits with A 7^ 0, as 
recognized in Fig. [TU1 fab 

In usual band structures, we have shown that dominant 
contribution arises from the Fermi surface term ail, and 
the relation given by eq. (|44p holds well. On the other 
hand, in the exceptional case such as in Ta, the anoma- 
lous 7-dependence of SHC can be explained by analyzing 



zllb 



as follows: By dropping the current operators in 
eq. (fTT)|) for simplicity, ^x y Ih is given by 
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We can approximate as follows for 7 <C \E l k \, \E k n \: 
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(47) 



Substituting above equation into eq. (|46p . we obtain the 
following relation for small 7: 
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k,l>rn 
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(48) 
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where the step function of the second term in above equa- 
tion is introduced to skip the fc-summation in the case of 
7>\E l k l\El n \. 

Figure [TO] (b) is a schematic band structure around the 
accidental degenerate point slightly away from the Fermi 
level. In this figure, A represents the band-splitting mea- 
sured from the Fermi level and E* represents the eigenen- 
ergy measured from the Fermi level: E* — E l k , ~ E^l . 

The first term of eq. (|48|) is finite only in region (B) 
where E l k > and E^ < 0, and its sign is positive The 
sign of the second term is negative in region (B), whereas 
it is positive in regions (A) and (C). When 7 — ► 0, o z J y Ih 
is given only by the first term in eq. ([48]) since the sec- 
ond term vanishes. When E* ~ 7, a large contribution 
to cr xy ° comes from the second term in eq. (|48p from the 
regime (C) in Fig. [TO] (b). Since the second term can be 
as large as the first term, cr%U takes the sizable peak at 
finite 7 in the presence of almost degenerate anticross- 
ing points near the Fermi level. The SHC reaches the 
maximum value at 7 <~ E* . As a result, the anomalous 
behavior of SHC in NRL-TB model for Ta originates from 
the anticrossing points as shown in Fig. [TO] (a). From this 
figure, we see that EL, is ~ 0.01. This fact is consistent 
with the peak of SHC a% y around 7 = 0.02 in Fig. (a). 

In the present model with A=0, there is another ac- 
cidental degenerate point at (k x , k y , k z )=(n, 0.237T, 7r/2). 
The band structure obtained for NRL-TB model in Ta 
with A=0 and A=0.023 around this point is shown in Fig. 
[TT1 From this figure, when /i is lowered to A/i = —0.04, 
we see that the Fermi level begins to lie inside the gap. 
This fact causes a sharp peak of <J xy Ih at Afj, = —0.04 in 
Fig. [8](b). 



0.05 
0.0 



5 



-0. 

















\ftx,kz)=(7t,7t/2)" 



nil 



FIG. 11: Anticrossing bands near the Fermi level with A = 
0,0.023 in Ta: The ^-dependence of E l k for (k x ,k z ) = 
(tt,tv/2) is shown. In the present model with A = 0, we find 
an accidental degenerate point at (k x , k y , k z )=(n, 0.237T, 7r/2), 
which is different point from that shown in Fig. \W\ (b). 
The corresponding minimum bandsplitting induced by SOI 
is - 0.015. 

Here, we have shown that the presence of almost degen- 
erate anticrossing points in Ta gives rise to the anoma- 
lous 7-dependence in the low resistive regime based on 
the NRL-TB model. In this exceptional situation, <J z J y Ih 
plays an significant role. Except for this special case, 
however, SHC is mainly given by the Fermi surface term 



a xy . We may have to confirm this anomalous behavior 
in Ta by checking the accuracy of the NRL-TB model in 
detail. 



V. DISCUSSIONS 

A. EFFECTIVE MAGNETIC FLUX 

In the previous section, we have discussed the SHC 
based on multiorbital tight-binding model using the 
Green function method. In this section, we give an intu- 
itive reason why huge SHC appears in the present mul- 
tiorbital model based on the double layer bec model in 
Fig. [TO] which is a simplified version of the bec struc- 
ture model. Here, we consider only d xy -, d x i_ y i- and s- 
orbitals considering the fact that d xy - and <i,.2_ a 2-orbitals 
give the dominant contribution to the SHC in various 
transition metals as explained in section IIVI In Fig. 
I12[ ±t represents the hopping integrals between near- 
est neighbor d^-orbital and s-orbital, and ±£' is for the 
next nearest neighbor d x 2 _ y i -orbital and s-orbital. Note 
that both hopping integrals change their signs by rota- 
tion by 7r/2. to represents the hopping integral between 
s-orbitals. 



s-orbital 



SO interaction 

\l s, : s=HI2 




s-orbital 












2 # ..;.;;;:/j£' 






3 5 








» 





AB phase factor 

±1 - e 
= ±0 o /4 



FIG. 12: Effective AB phase in double layer bec model for j- 
spin electron. This is the origin of the huge Hall conductivities 
in various transition metals. 



First, we explain that electron can transfer from 
dj2_ y 2-orbital state to d^y-orbital state and vice versa 
by using SOI for f-spin electron HXl z /2, which plays a 
significant role in the large SHE and OHE in transition 
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metals [Sj. \x 2 — y 2 ) state is transferred to \xy) state by 
rotating the wave function around z-axis by 7r/4. Since 
l z is a generator of the rotaion operator about z-axis 
R z {9) = e~ l91 * , the relation between l z and R z (8) is given 
by 



fl.(±i). 
for l a = ±2 basis. Therefore, 

Z> 2 - y 2 ) = 2 t R z m\x 2 - y 2 ) = 2i\xy), 



l z \xy) = -UR z {\)\xy) 



-2i\x 2 



y 



(49) 

(50) 
(51) 



As a result, we obtain the off-diagonal matrix element of 
SOI for f-spin electron as follows: 

(xy\hXl z /2\x 2 - y 2 ) = -(x 2 - y 2 \h\l z /2\xy) = ih\. 

(52) 

Figure [12] shows the most important process (interor- 
bital hopping process) for SHE in real space. By con- 
sidering the signs of interorbital hopping integrals (±£ 
and ±t') and matrix elements of SOI in eq. (|52|) . we 
can verify that a clockwise (anti-clockwise) movement of 
a t-spin electron along the path "0^3^1 or 2^0", for 
example, causes the factor +«(— i). This factor can be 
interpreted as the Aharonov-Bohm phase factor e 27 ™*/^ 
[4>o = hc/\e\], where (f> represents the effective magnetic 
flux cj> — § Adr — ±<fio/4. Since the effective magnetic 
flux for J,-spin electron is opposite in sign, electrons with 
different spins move to opposite direction. Therefore, 
the effective magnetic flux gives rise to the SHC of or- 
der 0(A). This mechanism will be realized in various 
multiorbital transition metals @, § ■ 



B. Coherent-Incoherent Crossover of Intrinsic Hall 
Conductivities 

In section ITVl we studied the 7-dependence of the SHC 
in Ta and W numerically. Therein, we have verified that a 
typical crossover behavior of the intrinsic SHC at 7 ~ A 
is realized in many transition metals unless almost de- 
generate anticrossing points exist slight away from the 
Fermi level. These crossover behavior of a xy was shown 
in ref. [H|,[22j]. Here, we discuss the crossover behavior of 
intrinsic Hall conductivities analytically by dividing into 
3 regimes with respect to 7: 7 « A, A « 7 < If , and 
W -C 7, where A and W represent the band-splitting 
near the Fermi level and the bandwidth, respectively. 
The first regime corresponds to the low resistive regime, 
and the second to the high resistive regime. Here, we dis- 
cuss the regime W <gC 7 only briefly since the Ioffe-Regel 
condition 7 ~ W is violated. 

Now, we analyze eqs. ([To]) . (fT5j) and (fT!))) to obtain 
the 7-dependence of the Fermi surface term a xy and the 
Fermi sea term a xy . The 7-dependence of a xy is esti- 
mated by analyzing eq. (Tj"rJ]) as follows for 7 <C W, by 



TABLE III: 7-dependence of the Fermi surface term a^ y , 
Fermi sea terms <T xy a ,crl.y , and the longitudinal conductiv- 
ity o xx . 





7 < A (low p) 


A < 7 < W (high p) 


W < 7 


®xy 
Ila 

® xy 
lib 
xy 


7 U 
7° 
7° 


I- 1 

7° 
7° 


7" 3 
7" 1 

7- 1 


xx 


7" 1 


7- 1 


I-' 2 



considering the following relationship: 

1 1 . liK-Ek) 



Im 



(E l - i~f)(E m + ij) 



((4) 2 +7 2 )((£D 2 +7 2 ) 
(53) 

* -y(E?-E l k )(6(E l k )+6(EZ)) 
7 {E™-E k y +1 * 



(54) 



After fc-summation, eq. (|54|) is proportional to 7 for 
7 <C A, and it is proportional to 7~ 2 for A <C 7 <C 
W . On the other hand, in the regime W <§; 7, we can 
estimate eq. (|53|) as 



E 



7(4 - K 



7 



((4) 2 +7 2 )((£r) 2 +7 2 ) 1 4 



E(4 



EV 



(55) 



In a similar way, the 7-dependence of the Fermi sea terms 
&x y a and al y h can be estimated from eqs. (fT8| and (fl"9j) . 
respectively. The longitudinal conductivity a xx is given 
by 



2ttA ^ 

k 



Tr I J™ G J~ G 



[ 

\{JCG R JCG R + (R^A)} 



(56) 



&XX 

a. 



In Table HTT1 the 7-dependences of a xy , a xy a , <J xy b , and 

are shown. In the metallic systems, the relations 

xy ~ a xy (Fermi surface term) and |cr:L| ^> \<r xy \ hold 

well for a wide range of 7 since the Fermi sea terms cr£f, a 

x y 

and a xy b almost cancel each other j22[. Therefore, we 
discuss the 7-dependence of the Fermi surface term in 
detail. From Table IIII1 the Fermi surface term a xy is in- 
dependent of 7 in the low resistive regime, whereas cr xy 
decreases approximately in proportion to 7~ 2 in the high 
resistive regime. On the other hand, the longitudinal con- 
ductivity a xx decreases approximately in proportion to 
7 _1 in both low and high resistive regime. Therefore, the 
coherent-incoherent crossover behavior of intrinsic Hall 
conductivities at 7 ~ A are reproduced by considering 
the Fermi surface term a xy correctly, as reported in refs. 

MM- 



a 



oc a 



a xy (X a 



XX 

2 

XX 



for 7 < A, (57) 
for A < 7 < W. (58) 
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In the case of W <C 7, the relation a xy oc <j x ^ (oc 7 -3 ) 
holds. However, this relation is not reliable since the 
Ioffe-Regel condition (W/j ~ Eft < 1) is violated in 
this regime [52| . 

Finally, we comment on the a xy b term, which is called 
the Berry curvature term. In electron gas models, the re- 



lation a xy — ,, xy 
M, IH. However, 



ui 1 } holds for 7 = +0 since a*. y + a 1 ^ = 

7xy 7^ a xy b f° r finite 7, and the 
crossover behavior cannot be explained by analyzing 
a xy b > which is shown in Tabic ITTT1 In conlusion, the re- 

olA > \all\ hold well in the 



a L xy and 



lationship a xy 

real metallic systems, and the correct crossover behavior 
given by eqs. ([57]) and (|58|) are reproduced by the Fermi 
surface term ai,,. 

y 

Now, we comment on the Ioffe-Regel limit in transition 
metals: The Ioffe-Regel limit l/a ~ kpl ~ 1 is approxi- 
mately estimated as Ep/j ~ 1, where /, a, kp, Ep, and 7 
represents an elastic mean- free path, a lattice constant, a 
Fermi wave number, a Fermi energy, and a quasiparticle 
damping rate, respectively. From the band structure for 
Ta and Pt in the present model, we verified that Ep is 
~ 1 in Ta and ~ 0.5 in Pt. Therefore, the Ioffe-Regel 
limit lays around 7 ~ 1 in Ta and 7 ~ 0.5 in Pt, respec- 
tively. Since the localization effect of electron is not taken 
into account, the present calculation will be inadequate 
for Ep/j > 1. 



VI. SUMMARY OF THE PRESENT STUDY 

In this paper, we studied the intrinsic SHE and OHE 
in various 4d-transition metals (Nb, Mo, Tc, Ru, Rh, Pd, 
and Ag) and 5d-transition metals (Ta, W, Re, Os, Ir, Pt 
and Au) based on a multiorbital tight-binding model. We 
have derived the general expressions for the intrinsic SHC 
and OHC in the presence of overlap integrals given by eq. 
©. We found that the huge SHCs in Pt (5c! 9 6s 1 ) and Pd 
(4d 10 5s°) are positive, whereas the SHCs in Ta (5d 3 6s 2 ) 
and W (5d 4 6s 2 ) take large negative values. We also found 
that the SHC changes smoothly with the electron number 
n = n s + rid, regardless of the changes of the crystal 
structure. Among the 4c?- and 5d-transition metals, the 
magnitude of SHC in Pt shows the largest value in the low 
resistive regime. However, the magnitude of SHC in Ta 
and W exceeds that in Pt in the high resistive regime. 
Therefore, large negative values of SHCs in Ta and W 
will be observed even in the high resistive samlpes. In 
this paper, we also calculated the SHC for n = 7, 8 (hep 
structure). 

We also showed that the CVC due to the local impu- 
rity potential can be safely neglected in calculating the 
SHC and OHC in the present model. The obtained SHCs 
in various transition metals are sensitive to the changes 
of the chemical potential /i, which reflect the multiband 
structure around the Fermi level /i. This suggests that 
the intrinsic SHC can be controlled by composing alloys. 
As for 7-dependences of SHC and OHC, we obtained the 
coherent-incoherent crossover behaviors in many transi- 



tion metals by calculating both Fermi surface term and 
Fermi sea term on the same footing: a xy is independent 
of 7 in the low resistive regime where 7 -C A, whereas a z xy 
decreases approximately in proportion to p~ 2 in the high 
resistive regime. The physical meaning of the crossover 
behavior can be explained as follows: In the low resis- 
tive rigime, SHC is proportional to the lifetime of the 
interband excitation h/A since it is caused by the inter- 
band particle-hole excitation induced by the electric field 
0, 0, 0, EH HE H2] • However, in the high resistive regime, 
SHC decreases with 7 since the quasiparticle lifetime fi/7 
becomes shorter than h/A (l5l. I22I]. 

Here, we comment on the effect of Coulomb interaction 
on the SHC and OHC. In the microscopic Fermi liquid 
theory, the Coulomb interaction is renormalized to the 
self-energy correction and the CVC. The renormalization 

factor due to the self-energy, z = (l — ^ ex- 

actly cancels in the final formulas of the SHC and OHC, 
eqs. JIBJ), (HJ) and |Q11). As shown in ref. [l5j, the CVC 
due to the Coulomb interaction does not cause the skew 
scattering . Therefore, 7-dependences of intrinsic SHC 
and OHC are unchanged by the CVC due to Coulomb 
interactions. However, it is well-known that the CVC 
causes various anomalous transport phenomena in the 
vicinity of the mag netic quantum critical points (QCP) 
[53L IB"! . [55L IBH [Hzf . In the same way, prominent CVC 
near the magnetic QCP may cause novel temperature 
dependence of the SHC and OHC. This is an important 
future problem. 

Owing to the present study of SHE and OHE in vari- 
ous transition metals, it has been revealed that the huge 
SHE and OHE are ubiquitous in multiorbital d-electron 
systems. In i ]V A[ we have discussed that the origin of 
these huge SHE is the "effective AB phase" induced by 
the atomic SOI with the aid of interorbital hopping in- 
tegrals @, [H ■ The present study strongly suggests that 
"giant SHE and OHE " will be seen ubiquitously in mul- 
tiorbital /-electron systems with atomic orbital degrees 
of freedom. These facts will enable us to constract effi- 
cient spintronics or orbitronics devices made of transition 
metals. Furthermore, in /-electron systems, a larger SHE 
may appear compared to that of d-electron systems since 
the angular momentum of atomic orbital is larger and 
the band splitting near the Fermi surface is smaller. 

In the presence of anticrossing bands, huge SHC can be 
realized when the Fermi level lies inside the gap induced 



by SOI [3QJ, |31|. In this case, a. 



, „ - and a xy takes a 
large (and almost quantized) value. In the present study, 
however, we could not find any elemental metals in which 
the Dirac cone type band structure has a dominant con- 
tribution to SHE in the low resistive regime. The large 
SHCs in transition metals are mainly given by the Fermi 
surface term, <r xy . Therefore, the existence of Dirac point 
is not a necessary condition for large SHE: As shown in 
Fig. [51 the band structure where the band-splitting near 
the Fermi level is small is significant for large SHE. How- 
ever, only in the case of Ta, almost degenerate anticross- 
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ing points that exists slightly away from the Fermi level 
gives rise to an anomalous 7-depcndcncc in the low re- 
sistive regime. The anomalous 7-dependence of SHC in 
Ta may be realized. 

Finally, we discuss the quantitative accuracy of the 
obtained numerical results, which depends on the accu- 
racy of the band structure of the model near the Fermi 
level. SHC and OHC depend on the multiband structure 
near the Fermi level with small interband splitting A, 
and they are proportional to A -1 according to eq. §TE§. 
According to ref. the possible error in the NRL-TB 
is about 0.002-0.004 Ry, which is much smaller than A 
(A ~ 0.035Ry in Pt) that gives the minimum energy scale 
in the intrinsic Hall effect. Therefore, it is expected that 
the NRL-TB model is accurate enough to derive qualita- 
tively reliable results of SHC and OHC. Thus, the overall 
n-dependence of the SHC in Fig. [5] will be reliable. 
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APPENDIX A: DERIVATION OF EQ. (|38|) 

Here, we derive eq. (|38|) . n(q) is given by Fourier 
transform of electron number density n(r) as follows: 

n(q) = J drn(r)e~ lqr 

dr^(r)^(r)e- zqr 



EE 

h,h' q,/3 



where i/j(r) represents the electron field operator and this 
operator can be expanded with the atomic wave function 
4>ki{f) as follows: 



(A2) 



We used this relation to transfer from the second row to 
the third row in eq. (|Alj) . 

From Bloch's theorem, the atomic wave function can 
be rewritten as 



1 



ik-r 



u ka {r)e lk ' 



(A3) 



Uk 



a (r) =J2e iHRi - r) Mr - Ri). 



(A4) 



Then, it is straight foward to show that the equation in 
square bracket in eq. (|Aip is rewritten as 



S(k + q-k') / dru* ka (r)u k >p(r) 

J unit 

Therefore, eq. f|Al|) is rcwriten as 



(A5) 



E 

k,a,p 



dru* ha {r)u k+q ,p{r) 



_J unit 



C la C k+q,0- (A6) 



Here, we expand [ ■ ] in eq. (|A6|) in powers of q/2 (q 
\q\) as follows: 



dri 



'k+q/2-q+2 



(r) (A7) 



= / dru L u k + U^ + o(l), (A8) 

J unit z vz/ 

where k is given by k = k + q/2, and §j4 Q( 3 is given by 
/ ^ 

J unit ■ ■ 



■■(Ri- r)e- fk ^- r ^* a (r - R t )e fk ^~ r Up(r - Rj) 

-| • (R j - r)e-^ R '- r U* a (r - i^e^-^r - Rj) 

(A9) 



First, the first term on the right hand side in eq. (|A8[) 
can be calculated as follows: 



dr Xy^^aCr - R i y k{R >~ r) Mr - fl,-) 



= o-Uk). 



(A10) 



Next, we show that \A a p given by eq. (|A9|) vanishes 
by rewritting eq. (|A9[) as follows: 



(Al) J2 S kk- 



q ( d d 



2 \dk dk 

k' 

dr^2e- fk ^~ r m{r - R t )e lk '^ R ^ r Up(r - Rj) 

(All) 



q ' to-l-^o-J)=o. 



2 a \dk ati dk Q/3 



Therefore, the final result for n(q) is given by 

«(<?) = E ^)^_ q/2 a c k+q/2 ^. (A12) 

k,a,P 

which is exact up to O(q). 
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APPENDIX B: DEFINITIONS OF SPIN AND 
ORBITAL CURRENT OPERATORS 

In the present study, we assume that the spin and or- 
bital current operators are given by eqs. (fTTj) and (JT3J) 
according to literatures 0, |47]]. Here, show the validity 
for these definitions in a microscopic way. Since the spin 
and the orbital operators are not conserved in the present 
model with SOI, we can not define spin and orbital cur- 
rent operators from the continuity equations. However, 
it is possible to make a natural definition for each current 
operator, as follows. 

First, we consider the spin current operator. Since the 
SOI in the present model is local, we can virtually apply 
a magnetic field (vector potential) to "["-spin and |-spin 
electrons separately. Here, we denote the vector potential 
for j-spin and i-spin as Af and Aj,, respectively. By 
considering a transformation k a — > k a — eA CT , where a 
is a spin index, the ^-component of current operators for 
|-spin and |-spin electrons are given by 



the real space representation as 



tC 



V J C , 



(B4) 



Jim.jrri ~ e ( r i r j) ' tim.jm' { c i m c jm' c j m i c im)j 

(B5) 

where i is the position of ith lattice point, m represents 
the eigenvalue of l z , and ti m j m i represents the hopping 
integral between \m) state at ith site and |m') state at 
jth site, respectively. Then, the natural definition of the 
orbital current operator will be given by 



(J )i 



1 



(m + m')(J ) im ,j 



2(-e) 



(B6) 



in the basis of l z = 2,1,- • -,—2. By perfoming Fourier 
transforms of eq. (|B6j) . we obtain the following expres- 
sion for the orbital current operator in the present model: 



J[x — 



dH l 


' v x 







v 





dH | 


'0 





dA lx ~ e \ 


v° 


v x 



(Bl) 
(B2) 



Therefore, the natural definition of spin current operator 
in the present model is given by 



s 1 



(B3) 



This expression is equivalent to eq. (fTTj) . We note that 
above discussion can not be applied to systems in the 
presence of non-local SOI, such as Rahsba type SOI. 

Next, we consider the orbital current operator. The 
charge current operator J x = J-t x + J\ x is expressed in 



(J°(k)) r 



1 



-(m + m')(J c (fc)) 



(B7) 



2(-e)' 

In general basis, the above equation can be rewritten as 



{J C .l z } /2(-e). 



(B8) 



This expression is equivalent to eq. (|12p . In summary, 
we have introduced a natural definition of the spin and 
orbital current operator, and shown that they are equiv- 
alent to eqs. (fTTj) and (fT2"|) , respectively. 

In the same manner, if we define the spin cur- 
rent operator (in the presence of the intersite SOI) as 

{J ^)ima,jm' a' = 2( — e) \* ® )\^ )im(T,jm' u' i Can de- 
rive the spin current J = {J c , s z } /2(— e) immediately. 
This is another microscopic derivation of the spin current 
operator in eq. (fTTj . 
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